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Влияние кинематических параметров движения локтя на электромиографический сигнал 
двуглавой мышцы плеча” 


Ф. Бонилья, Е. А. Лукьянов, А. В. Литвин, Д. А. Деплов 


Рассмотрено влияние кинематических параметров движения локтя (угловое перемещение, скорость и ускорение) 
на поверхностный электромиографический сигнал двуглавой мышцы плеча при выполнении движений локтя с раз- 
личной нагрузкой. Установлены зависимости частотных и амплитудных параметров ЭМГ-сигналов от параметров 
движения локтя при подъеме и опускании груза. Показана структура экспериментального стенда. Представлена 
методика записи ЭМ/-сигналов с двуглавой мышцы плеча, а также методика обработки сигналов. Для формирова- 
ния количественных значений показателей анализа сигналов использовались методы спектрального и статистиче- 
ского анализа. Статистический анализ во временной области включал в себя определение следующих параметров: 
дисперсия амплитудных значений ЭМГ, среднее и среднеквадратическое значение абсолютных величин амплитуд 
ЭМГ, частота пересечения нулевой линии ЭМ/-сигналом. Эти параметры, а также параметр, характеризующий 
форму сигнала, определялись на основе массива измеренных значений сигнала, формируемого в течение (0.3 се- 
кунды текущего интервала времени («скользящий» интервал). Оценена величина развиваемого мышечного усилия, 
найдены коэффициенты косинусного преобразования Фурье. Построены гистограммы распределений биопотенци- 
алов ЭМГ, В частотной области выполнялся спектральный анализ ЭМ/-сигналов методом быстрого преобразования 
Фурье. Определялись общая мощность спектра, средняя частота, медиана, частота спектральной составляющей, 
имеющей наибольшую амплитуду. Регистрация и анализ биопотенциалов ЭМГ выполнялись средствами Майа. 
Выявлены информационные признаки, которые могут быть использованы для синтеза интеллектуальной системы 
управления на основе нейронных сетей. 
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Введение. Скелетные мышцы являются совокупностью двигательных единиц (ДЕ). При стимуляции 
нейронным сигналом каждая ДЕ сокращается и генерирует электромиографический (ЭМГ) сигнал, ко- 
торый представляет собой сумму потенциалов действия всех вовлеченных в процесс клеток [1]. По- 
верхностная ЭМГ часто используется для оценки относительного уровня мышечной активности во время 
движений [2]. Известно, что слабые волевые усилия заставляют ДЕ возбуждаться с частотой примерно 
5—15 Гц, а при возрастании усилий частота увеличивается до 25-50 Гц с формированием интерферен- 
ционного типа ЭМГ-сигнала. По мере утомления мышцы наблюдается снижение доли высокочастотных 
составляющих и увеличение амплитуды ЭМГ [3]. Сигнал ЭМГ, записанный с использованием поверх- 
ностных электродов, является сложным, включает в себя интерферирующие составляющие, поэтому 
его анализ затруднен [4]. 

ЭМГ-сигнал характеризует уровень активности мышцы и широко используется в медицине для 
диагностики нейромышечных заболеваний. Достаточно активно ведутся работы по использованию ЭМГ- 
сигналов для био- или нейроуправления протезами [5-7]. Очевидно, для создания эффективных про- 
тезов с биоэлектрическим управлением необходим качественный и стабильный электромиографический 
сигнал. Также важно использовать эффективные усилители биопотенциалов, определить точки уста- 
новки электродов [7]. В статье приведены методика и итоги экспериментальных исследований электро- 
миографических сигналов, результаты их анализа, позволяющие формировать набор признаков для 
распознавания выполняемого движения и его параметров. 

Методика проведения исследований. Исследования проводились при помощи специальной элек- 
тромеханической системы типа экзоскелет, которая ограничивала подвижность локтевого и плечевого 
суставов (рис. 1). В исследованиях принимали участие пять испытуемых в возрасте от 24 до 35 лет, 
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примерно равные по массе тела и физическому развитию, не имеющие заболеваний и повреждений 
верхних конечностей. 

Канал записи электромиографического сигнала включал накожные электроды, усилитель, ана- 
лого-цифровой преобразователь (АЦП), программные фильтры: режекторный 4-го порядка (для удале- 
ния сетевой помехи) и полосовой с частотами среза 10 и 250 Гц. 

Сигнал ЭМГ записывался системой сбора данных МТ У$В 6212 с помощью электродов Ад/АдС/ 
(диаметр 5 мм). Частота оцифровки — 1000 Гц. Место установки электродов предварительно обраба- 
тывалось в соответствии с рекомендациями [4]. 

Электроды крепились на поверхность бицепса по линии мышечных волокон на расстоянии 2 см 
друг от друга. Они фиксировались с помощью манжеты, обеспечивающей их надежный контакт с кожей 
предплечья. Это исключало перемещение электродов в процессе эксперимента (рис. 1). 





то 





а) 6) 
Рис. 1. Электромеханическая система для записи ЭМГ-сигнала бицепса: общий вид экзоскелета (а); 
кинематическая схема (6); 9 — угол сгибания локтя 


Цель исследования — установление зависимостей характеристик ЭМГ-сигнала бицепса от пара- 
метров движения локтя. Кроме того, выявлялись характерные параметры, позволяющие определять 
фазы движения локтя для управления внешними устройствами. 

Значения углового перемещения плеча (см. рис. 1) измерялись с помощью гироскопа (датчик 
[3 642000). Данные (как и ЭМГ-сигнал) записывались в файл на компьютере. Для этого использовалась 
система сбора данных МТ Ц$В-6212. Испытуемый изменял угловое положение плеча в диапазоне 0-90-— 
0 градусов при нагрузках 5, 7, 15 кг. 

Анализ данных. Полученные экспериментальные данные исследовались методами статистического, 
спектрального и вейвлет-анализа. Использовалось также авторегрессионное моделирование [1]. Для 
исследования электромиографических сигналов во временной области применены скользящие окна, 
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размер которых определялся условиями эксперимента. После записи сигнала определялись его диспер- 
сия, среднее абсолютных значений (/77 АИ), среднеквадратическое значение ( АЛ/5), частота пересечения 
нуля (ГС), параметр, характеризующий «форму» ЭМГ-сигнала (Иеп). 

Развиваемое мышечное усилие оценивалось с учетом параметра /одреесЕ а также коэффици- 
ентов косинусного преобразования Фурье [1, 8]. 

Дисперсия является усредненной характеристикой вариабельности ЭМГ-сигнала относительно 
его среднего значения. Скользящее среднеквадратическое значение — индикатор средней мощности 
ЭМГ-сигнала как функции времени. Выбирался такой размер окна, чтобы на основе анализа имеющихся 
в окне данных можно было принимать решения о выполняемом движении. Экспериментально опреде- 
ленная величина окна составила 0,2-0,3 секунды. 

Значение показателя ГС определяется количеством пересечений ЭМГ-сигналом нулевой линии 
в пределах окна. Этот показатель является простой мерой, связанной с частотой и сложностью сигнала. 
Для устранения влияния низкоамплитудных шумов при обработке сигнала вводился порог, равный 10-— 
15 мВ. 

Параметр И/еп (И/ауеюгт [епдё”), характеризующий форму сигнала, является мерой его 
«сложности» и определялся на основе отсчетов в текущем окне [8] следующим образом: 





№ 
И/ауеет =У`|^х, |, где Ахх = Хк — Хи. (1) 
К=1 


Развиваемое мышечное усилие оценивалось посредством параметра /одреесЁ по следующей 
зависимости: 


(им У юар 
ЮдреесЕ =е “м (2) 


Коэффициенты авторегрессионных (АР) моделей сигналов ЭМГ, содержащие информацию об 
интенсивности сокращения мышцы [8], определялись из выражения: 


р 


где а/ — коэффициенты АР-модели; р — порядок АР-модели; е— остаточный шум. 
Коэффициенты косинусного преобразования Фурье [8] определялись из авторегрессионных мо- 
делей: 


/—1 
с =-а; а, =а, - > Па ы (4) 


где а; — /й коэффициент АР-модели, с — /-й коэффициент косинусного преобразования Фурье. 

Определялись также стандартные амплитудные параметры, такие как среднее значение пиков, 
максимальное и минимальное значение пиков, площадь ЭМГ-сигнала. 

Для спектрального анализа ЭМГ-сигналов использовался метод быстрого преобразования 
Фурье. Определялись общая мощность спектра, средняя частота, медиана, частота спектральной со- 
ставляющей, имеющей наибольшую амплитуду [4, 9]. 
Результаты исследований. Проведенная работа позволила определить зависимость параметров 
ЭМГ-сигнала биопотенциалов бицепса от углового перемещения и скорости плеча (табл., рис. 2-3). В 
таблице верхнее значение параметра соответствует подъему плеча, нижнее — опусканию. 

В таблице приведены результаты, рассчитанные для ЭМГ-сигнала, ограниченного окном в 300 
мс в моменты времени, соответствующие максимальным скоростям движения плеча. 

Из таблицы видно, что при подъеме/опускании значимо отличаются параметры /одреесёи ГС. 
Так, при нагрузке 5 кг и времени выполнения движения 6 с для параметра /одреесЁ отличие составило 
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32 %, для параметра 2{С— 42 %. При нагрузке 15 кг и времени перемещения 6 секунд соответствующие 
отличия для /одреесЁи ГС— 7,9 % и 18,9%. 

При нагрузке 5 кг увеличение времени перемещения до 8 с привело к снижению различий 
ЮЮдрЕесЕ — 6,3 %, т. е. параметр при опускании был несколько больше сигнала при подъеме. Это 
объясняется физиологией изменения ЭМГ-сигнала при длительном мышечном напряжении. При этих 
же условиях параметр 2С изменяется при подъеме/опускании на 38 %. Увеличение нагрузки до 15 кг 
(время 8 с) при сгибании и разгибании локтя обеспечивает изменение параметра /одреесЁна 25 %, ГС 
на 8 %. 


Таблица 1 


Зависимость параметров ЭМГ-сигнала биопотенциала бицепса от силовой нагрузки 
и кинематических параметров движения локтя 


град/с 

67,09 0, Е 0,043 Ч С 6,786 
Я ИБьИ -ШЕи-ЖЕЕ- 

59. 50 о, 0514 о, 034 о, - 5,161 
7” ыы ом | 0 | ив 

56 92 о, 0489 0, 032 0, о 4,529 
о о | ом мы 


\Меп, мс 


Вес, | Время, с| Длина окна, \/тах, |од)еес Р 
с Пььы 

58,54 0,2233 0,163 0,29 28,11 
авиа 


53,29 0,2081 0,132 0, 36,88 

ИЕ ИЕ Е 

23, 59 б, 2965 б, 156 б, > 45,04 

т 5 о | о [0% 

Из рис. 2, 3 видно, что наибольшие величины ЭМГ-сигнала соответствуют максимальному раз- 
виваемому усилию. 

Во время выполнения экспериментов отмечено, что амплитудные значения ЭМГ-сигнала при 
сгибании и разгибании локтя неодинаковы (при сгибании больше) и зависят от нагрузки. При нагрузке 
15 кг сигнал в 5,9 раза больше, чем при нагрузке 5 кг. Время, за которое совершались движения плеча 
при сгибании и разгибании, определялось по осциллограмме и составляло примерно 4 си 3,9 с 
(рис. 3 в, Г). 

Анализ в частотной области показал, что увеличение нагрузки сопровождается увеличением 
частоты разряда каждой двигательной единицы [1]. Поднимание груза 5 кг вызывало появление мак- 
симальных частот в интервале 25-50 Гц. При нагрузке 15 кг интервал частот смещался в область 
40-65 Гц (см. рис. 3, а, 6). Максимальная величина мощности спектральной плотности отличается в 
20 раз. Аналогичные результаты получены и для спектрограмм (см. рис. 3, в, г). Это свидетельствует о 


том, что при возрастании нагрузки формируется интерференционный тип ЭМГ, а пространственно-вре- 
менное суммирование всех активных двигательных единиц вызывает увеличение сигнала ЭМГ [1]. 
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Анализ данных во временной области показал, что изменение параметров /одреесЁ ГСи Иеп 
зависят от направления движения. Эти параметры в основном уменьшаются в фазе разгибания локтя. 
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Рис. 2. Влияние углового перемещения и скорости локтевого сустава при поднимании и опускании груза 5 кг 
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в течение 7 с: угловое перемещение и скорость (а); ЭМГ-сигнал бицепса (6); 
спектр мощности сигнала ЭМГ (в); спектрограмма (г) 
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Рис. 3. Поднимание и опускание груза 15 кг в течение 7 с: угловое перемещение и скорость (а); 
ЭМГ-сигнал бицепса (6); спектр мощности ЭМГ-сигнала (в); спектрограмма (г) 


На рис. 4 показано, что можно построить гиперплоскость для разделения на два класса значе- 
ний параметров 2Си /одрееа ЭМГ-сигнала бицепса при движении локтя вверх и вниз. 
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Рис. 4. Скаттерграммы 2С(а) и юдреесЕ (6): значения 2Спри движении локтевого сустава от 09 до 859 отмечены 
кружками; значения /одреесЁ при движении локтевого сустава от 859 до 09 отмечены звездочками 


Заключение. Исследование влияния кинематических параметров движения локтя на поверхностный 
ЭМГ-сигнал двуглавой мышцы плеча было выполнено в соответствии с современными рекомендациями 
и требованиями. Полученные результаты сопоставимы с данными из информационных источников. Вы- 
явленные информационные признаки могут быть использованы для синтеза интеллектуальной системы 
управления на основе нейронных сетей [10], которая позволит распознавать ЭМГ-сигнал и определять 
разные фазы движения. На этой основе может быть организовано управление, например, исполнитель- 
ными электромеханическими системами в условиях неопределенности. 
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Технические науки 


ЕРРЕСТ ОЕ КТМЕМАТТС РАКАМЕТЕВ$ ОЕ ЕЕВОМ/ МОТТОМ ОМ ВТСЕР$ ЕСЕЕСТВОМУОСКАРНТС 
$ТСМАЕ^ 


Е. ВопШа, А. Е. ЦиКуапом, А. \. мт, О. А. Оер]оу 


ТЛе еПесЕ о! те еБои/ тойоп Ктетайс рагатеег$ (апди!аг о@р/асетепЕ уе/осйу, апа ассеегайоп) оп {Пе зигГасе 
е/есготуодгар!/с (ЕМС) хата! о! пе 6/сер$ о! {пе агт ипдЕег те еБои/ товоп ехесийоп ив а!ТегепеЕ /оад /5 апа/угед. 
Тле аерепдепсез о! {пе Гедиепсу апа атрШиае рагатеег$ о! ЕМС; 51дпа/5 оп {пе еБои/ тойоп рагатеег$ ипдЕег /оаа 
Гепа апа Гоиегпад аге е$аЬ/ред. Тре (е5Е Бепсй $гисёиге 15 пои. Тре (еспт/дие о! гесогата ЕМС; $/дпе/5 гот пе 
Б/сер$, а5 и/Е// аз ОР (Пе 5/дпа! ргосЕ$$та теосб аге ргезегед. Меосб о! а 5йса!/ апа 5ресёга! апа/у5/5 аге и5ЕС Гог 
те диапИйсайоп о! пе ЕМС. 51апа! уа/иез. п пе ите аота!, те %айв5йса! апа/у$/5 псиде$ те деегт/пайоп о! пе 
юЮ/Юи/па рагатеегз; ЕМС атр/иае иа/иез орег5/оп, ауегаде апа КМ5` иа/иез, гего сгоз$/па иа/ие, игауегогт епдоЕв. 
ТЛе5е рагатеЕг$ аге деегт/пед оп {Пе теазигед иа/ие$ а/ггау Ва$/5 ог {Пе х/дпа! депегаед ии! 0.3 $ес. ог пе сиггепё 
ите (51Ата’ 'егуа/). ТВе деие/ореа тизсЕ Гогсе иа/ие 15 езитаед, те сера! сое"сепЕ5 аге Гоипа. ЕМС' БГороеп- 
ва! азтриноп Риодгат5 аге сопугисЕед. [п пе едиепсу адотат, пе ЕМС, 5/дпа!$ зресга! апа/у$/5 /5 сагтед сиё Бу 
Еа5Ё роитег ТгапзГогт. Тве 5ресгит а! рои/ег, йе ауегаде Гедиепсу, те тед/ап, апа а ресёга/ сотропепЕ ий” 
те аот/папЕ гедиепсу аге аЕ®егттпеа. Тре ЕМС Б/ороёепба! гед/5гайоп апа апа/уз/5 аге регГогтед бу теап$ о! Майа. 
[пГогтайоп Геайигез и/Р/СЙ сап бе изед Гог {ве е//депЕ сопёго! сует 5уп е5/5 Базеа оп пеига! пе //огк$ аге Аепийесц. 
Кеуигога$: ра55/уе ехозке/еюп, еБои/ товоп рагатеЕег5, еесготуодгарйи/с 5/дпа/, $ ай5йса! апа/у5/5, 5ресёга!/ хата! 
апа/у$/5, Майа. 


“ ТВе гееагсй 15 допе ут Не #гате оЁ Не тдерепаепе В&р. 
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